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Abstract 



The band gaps and spectral shifts of CdS, CdSe, CdTe, A1P, GaP, GaAs, 
and InP semiconductor clusters are calculated from band structure calcula- 
tions using accurate local and non-local empirical pseudopotentials. The effect 
of spin-orbit coupling on the band structures is included in the calculations 
when they are important. The complete set of pseudopotential parameters 
and full computational details are reported for all these semiconductors. The 
calculated spectral shifts of zinc-blende and wurtzite CdS, wurtzite CdSe, 
and zinc-blende InP clusters are in good agreement with experiments over 
a range of cluster sizes. The effect of crystal structure on the band gaps is 
small in large clusters but becomes important in small clusters. In the ab- 
sence of experimental data, our calculations provide reasonable estimates of 
expected spectral shifts for the other clusters. These results demonstrate that 
the empirical pseudopotential method yields unique insights into the quantum 
confinement effects and is a powerful tool for calculating the spectral shifts of 
semiconductor clusters. 

PACS numbers: 71.35. +z, 36.20.Kd, 36.40. +d 
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I. INTRODUCTION 



Clusters are an embryonic form of matter whose microscopic study provides insights into 
the evolution of material properties from molecules and surfaces to solids []T|— . Further- 
more, clusters have been shown to exhibit exotic optical properties and reactivities quite 
different from those in molecules and solids For these reasons, theoretical studies on 

clusters are critical to the design and synthesis of advanced materials with desired optical, 
electronic, and chemical properties |?HT6[. Such studies are at the interface of the tradi- 



tional fields of quantum chemistry, solid state chemistry, and statistical mechanics fi7H20 



Hence, physicists, chemists, and material scientists are working individually and in teams 
to unearth the fundamental principles underlying the structure, dynamics, and reactivities 
of these clusters pT|-f2(j|]. However, theoretical studies on the spectroscopy of semiconductor 
clusters have lagged far behind [ p?fl . 

The interest on the spectroscopy of semiconductor clusters arose from the discovery 
by Louis Brus that such clusters can be synthesized in colloidal suspensions by controlled 
liquid phase precipitation reactions |28| . The radii of the clusters thus synthesized are 
usually in the 5-100 A range. Furthermore, the X-ray and transmission electron microscopy 
(TEM) experiments have shown that these clusters have the same lattice structures as the 
corresponding bulk materials QJ^,^,^,^]. Finally, it is possible to prepare macroscopic 
samples of these clusters either in powder form or in colloidal suspension form. Such a 
versatility in sample preparation has made it possible for experimentalists to investigate 
the physical and optical properties of these clusters in detail f|. Indeed, very sophisticated 
experiments are currently underway in various laboratories to understand the absorption 
spectra of these clusters as a function of the cluster size |5, 17-19]. What is remarkable 
about these clusters is that the electronic spectra of these clusters are not the same as that 
of the bulk, even though the clusters have the same structural properties as the crystals 
Q. The absorption spectra of clusters show relatively sharp resonances superimposed on 
a continuum. The first absorption peak in the spectrum corresponds to the threshold for 
the absorption of light by the semiconductor cluster. It corresponds to the exciton energy, 
the energy needed to excite an electron from the top of the valence to the bottom of the 
conduction band. Quantitative prediction of the shift of the exciton energy with cluster size 
had been an outstanding problem for a decade. Simple theoretical calculations based on an 
effective mass model (EMM) were not successful and ab initio quantum chemistry electronic 
structure calculations are impossible for these clusters consisting of thousands of electrons. 
Consequently, what was needed was a simple model that captured the essential physics of 
the problem at hand. 

Recently, we developed such a model and solved this problem J3"0"| . We used accurate 
pseudopotentials to carry out band structure calculations and obtained the electronic energy 
levels of these clusters. Our calculations of the threshold for the absorption of light by CdS 
clusters yielded results in excellent agreement with experiment over a range of cluster sizes 
3D|] . Furthermore, we predicted two new effects. First, we found that the exciton energies 



in small CdS clusters are sensitive to the crystal structure, even though such a sensitivity is 



absent in large clusters |30|]. Second, we found that the vertical Franck-Condon transition 



energies in indirect gap clusters exhibit an anomalous redshift in small clusters [30]. In sharp 



contrast, the extant theoretical models always yielded monotonic blueshift of the transition 



2 



energies with decreasing cluster size |30] 



The band structure model we employ here has also yielded a noteworthy insight into 
the nature of the electronic transitions in semiconductor clusters. The question is whether 
the electronic transitions observed in these clusters are related in any way to the bulk tran- 
sitions. Previous theoretical calculations implicitly assumed that the electronic states in 
semiconductor clusters are similar to those in molecules and hence they are qualitatively 
different from those found in bulk. In sharp contrast, our band structure model underlies 
the presumption that the electronic states in nanoscale clusters are similar to those in bulk. 
Consequently, we use the bulk language (direct vs. indirect) to classify the electronic tran- 
sitions in clusters. Recent experiments of Brus and co-workers [f29| on the luminescence of 
Si clusters support our assumption. Based on the analysis of the luminescence of Si clus- 
ters, Brus and co-workers concluded that the lowest energy exciton transition is essentially 
an indirect gap transition, eventhough the spectrum has blueshifted by almost 0.9 eV ||29| . 
This observation that bulk-like transitions are preserved even in clusters supports our band 
structure model. 

In this paper we extend our band structure calculations to a variety of semiconductor 
clusters for two reasons: First, we wish to investigate the reliability of our band structure 
model by repeating these calculations on different systems for which accurate experimental 
data are now available. Second, the accuracy of our calculations on wurtzite clusters could 
not be verified previously because of the absence of experimental data on those clusters at 
that time. Since highly accurate experimental data are now available for wurtzite CdS and 
CdSe clusters, we felt the need to repeat our band structure calculations on these clusters. 
In this paper we also give the results of our calculations on CdTe, A1P, GaP, GaAs, and InP 
clusters. While accurate and reliable experimental data on these clusters are not available at 
present, we hope that our calculations provide reasonable estimates of the expected spectral 
shifts function of cluster size. 

Experimental synthesis of semiconductor clusters is still a challenging problem at present. 
One major difficulty is establishing a suitable synthetic route that yields nanometer scale 
clusters. The other major difficulty is determining optimal control parameters that yield 
clusters with a narrow size distribution. While both these problems remain to be fully 
solved, significant progress is being made daily. We hope the present paper will serve as a 
useful guide to the experimentalists on the expected behavior of spectral shifts in a variety 
of technologically important nanoscale semiconductor clusters. 

This paper is organized as follows: Section II gives the empirical pseudopotential theory 
of band structure calculations, Sec. Ill presents the band structure model, Sec. IV gives 
some computational details, and Sec. V presents results of these calculations on a variety 
of semiconductor clusters. Finally, we summarize this paper in Sec. VI. 



II. EMPIRICAL PSEUDOPOTENTIAL METHOD 

The electronic structure calculations of the semiconductor clusters are carried out us- 
ing the empirical pseudopotential method (EPM) that has been previously used for the 
investigations of the optical properties of bulk semiconductor materials and clusters [BO 



This method consists of solving the Schrodinger equation using an empirically determined 
pseudopotential for the valence electron, 
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(i) 



V p (v,E) = V L (r) + J2^ME)fi(mi, (2) 
1=0 

where the first term V L is the purely local part, the second term gives the non-local (VJvl) 
part, and Hi is the projection operator for angular momentum states I. The local part of 
the pseudopotential is given by 

Vl(t) = £ [V S (G)S S (G) + iV A (G)S A (G)} exp(*G • r), (3) 
G 

where Vs and Va are the symmetric and anti-symmetric form factors, respectively, deter- 
mined by fitting them to the experimental optical data. Similarly, Ss and Sa are the 
symmetric and anti-symmetric structure factors, respectively, determined from the crystal 
structure. The function Ai(E) is the energy dependent well depth 

ME) =a l + friE^K^iK')] 1 / 2 - E°(K F )}, (4) 

where E°(K) = h 2 K 2 /2m and K F = (6vrV^c) 1/3 is the Fermi wave- vector with z the num- 
ber of valence electrons per unit cell and Q c the volume of the unit cell. fi(v) is conveniently 
taken to be the square well 

M -{I III < 5 > 

ail, Pi, an d Ri are the non-local parameters of the EPM theory to be determined from the 
experimental optical data. 



A. Non-local Pseudopotentials 

In many instances the local part of the pseudopotential (Vl) is sufficiently accurate 
to represent the gross features of the band structures correctly. However, the non-local 
pseudopotential is much more accurate over a broad range of energy scales. In particular, 
the calculated band widths and band dispersion are closer to the experimental values with the 
non-local calculations than with the local calculations. The non-local calculations are also 
much more satisfactory theoretically since they incorporate the correct angular momentum 
dependence nature of the pseudopotential experienced by the valence electrons. 

The local contribution to the pseudopotential is given by Eq. (^). Evaluation of these 
matrix elements of the Hamiltonian will be described in Sees. IIC and IID for zinc-blende 
and wurtzite structures, respectively. In the plane wave basis the matrix elements of the 
non-local pseudopotential are of the form 

Att 

V NL (K,K') = —J2A}(E)P l (cos(e KK ,))S\K-K')F?(K,K'), (6) 

a Li 



as discussed elsewhere [31 . 
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B. Spin-Orbit Coupling in Band Theory 



In light elements the electron spin (s) and orbital angular momentum (I) are both good 
quantum numbers, since the magnetic field generated by the orbiting electron is too weak to 
induce coupling with the electron spin. However, in heavier elements the nearly relativistic 
speed of the valence electron produces a sufficiently large magnetic field that I and s are 
coupled, giving rise to j = I + s as the good quantum number. The spin-orbit interaction 
thus couples electron dynamics in spin and ordinary spaces, thereby reducing the overall 



symmetry of the Hamiltonian. This relativistic effect is represented by the operator [ 32 1 



where a are the Pauli spin matrices, V p is the pseudopotential, p is the electron momentum, 
m is the true electron mass, and c is the speed of light. The matrix elements of the new 
Hamiltonian in the plane wave representation are given by |33| 

< k + G', s'|H|k + G, s > = fi|k + G\ 2 5 GG >6 SS/ 

+ S S (G - G')[Vs(\G - Gf)5 ss , - z\ s (G' x G) • a ss ,} 

+ iS A (G - G')[V A (\G - G'\ 2 )5 SS , - i\ A (G' x G) • (8) 

For a binary semiconductor consisting of two types of atoms (A ^ B), we define 

A 5 = ^(Ai + A 2 ), A A = |(A 1 -A 2 ), (9) 

where 

Ai = ^(G)B^(G'), A 2 = a^Bf d (G)Bf d (G'), (10) 

fi is the adjustable parameter chosen in order to obtain the splitting A of the valence band 
at T correctly, and a is the ratio of the contribution from atom A to the contribution from 
atom B at T. The B n i are defined as 

/•oo 

B nl (G)=C r'RUrMG^dr, (11) 

<j 

where ji are the spherical Bessel functions, C is determined by the condition 

hm ^ = 1, (12) 
G^O G 

and R n i are the radial parts of the outermost electron wave functions taken from the Herman- 



Skillman tables [0. For simplicity we always used a value of a = 1.0 and C = 1.0. 

The band structures of semiconductors have several common features: At T the HOMO 
is split by the spin-orbit coupling Specifically, counting the spin, for the zinc-blende 
structure the HOMO band is sixfold degenerate. The spin-orbit interaction splits this band 
into an upper fourfold degenerate r 8 component, and a lower twofold degenerate T 7 com- 
ponent. The conduction band is twofold degenerate |36fl . For the wurtzite structure, the 
degeneracy of the T band is already removed by the crystal field, leading to an upper twofold 
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degenerate band, and a lower fourfold degenerate band. The spin-orbit coupling then splits 
the lower band into two twofold degenerate bands ]36|]. For both systems, wurtzite and 
zinc-blende, the highest three valence bands are called heavy-mass, light-mass, and split-off 
band, respectively []36[| . Transitions between these three valence bands to the conduction 



band are traditionally labelled A, B, and C ||37|| . The atomic orbital approach is especially 
useful for describing valence bands near T. From this perspective, the HOMO bands arise 
from the valence p orbitals while the LUMO band arises from the s orbitals [SB]. In polar 



semiconductors such as CdS and CdSe, HOMO bands originate primarily from the anion, 
and LUMO bands from the cation. However, in non-polar semiconductors such as CdTe and 
GaAs, considerable mixing of the cation and anion orbitals will take place. Since our interest 
is in reproducing the splitting of the valence p bands near HOMO, our calculations consid- 
ered only the contributions from the outermost p orbitals. We do not consider spin-orbit 
effects on the innercore or the <i-electron states since these states are eliminated in EPM 
through the use of the pseudopotential. The s orbitals do not exhibit spin-orbit splitting 
since their orbital angular momentum is zero. 



C. Local Pseudopotentials of Zinc-blende Structure Crystals 

The local pseudopotential (Vl) experienced by a valence electron at a point r inside a 
crystal lattice is given by 



Vb(r) = X;«i(r-R-d i ) ) 



(13) 



where the summation is over all the basis atoms j at each lattice point R, and v j is the 
atomic pseudopotential due to the jth basis atom at the lattice site R. Fourier expansion 
of Vj yields, 



Nn 



£5>,(G)exp[zG-(r-R-d,)], 
a G Rj 



(14) 



where N are the number of unit cells and n a are the number of basis atoms per unit cell. 

The zinc-blende (or sphalarite) lattice consists of two interpenetrating fee lattices, dis- 
placed from each other along the body diagonal by a /4, a being the lattice constant of 
the unit cell. Consequently, we may view the zinc-blende structure as an fee lattice with 
two different A ^ B basis atoms per unit cell. For this case, n a = 2,Vi(G) ^ vi(G), and 
di = — d 2 = — ti = —(1, 1, l)/8. Explicitly summing the above equation now yields, 



Vt.(t) 



E 

G 



1 



Vj(G) exp(— iG ■ dj 



exp(«G • r) 



[l MG) exp(*G • ta) + v 2 {G) exp(-iG • ti)]} exp(zG • r) 

G 

J2 [V S (G)S S (G) + iV A (G)S A (G)} exp(*G • r), 
G 



(15) 

(16) 
(17) 



where 
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V S (G) = \ K(G) + ^(G)] , V A (G) = \ K(G) - ^(G)] 

are the symmetric (Vg) and anti-symmetric (Va) form factors and 

S S (G) =cos(G-ti), S A (G) = sin(G • ti) 



(18) 



(19) 



are the symmetric (Sg) and anti-symmetric (5a) structure factors, respectively. Defining 
the reciprocal lattice vectors G as 



2n 

G = — (G x , G y , G z ) 
a 



we obtain 

Ss(G) 



cos 



71 



— (G x + G y + G z ) 



S A (G) = sin 



7T 



— (G x + Gy + G^) 



(20) 



(21) 



If we now specialize to the case of A = B, we obtain the local pseudopotential for diamond 
structure: 



V L {v) = J2 V S (G)S S (G) exp(*G • r). 



(22) 



G 



D. Local Pseudopotentials of Wurtzite Structure Crystals 



The wurtzite crystals are made of two interpenetrating hexagonal close packed (hep) 
lattices. One hep lattice is entirely made of A type atoms and the other entirely of B 
type atoms (A ^ B) and these two lattices are displaced from each other by 2t 2 = |co 
along the c-axis. However, the hep lattice is not a Bravais lattice. The hep lattice consists 
of two interpenetrating simple hexagonal Bravais lattices, displaced from one another by 
2t 1 = (|ai, |a 2 , |a 3 ), where (a 1 ,a 2 ,a 3 ) are the direct lattice primitive translation vectors 
of the simple hexagonal Bravais lattice. Hence, the wurtzite structure is a network of four 
simple hexagonal lattices, with four atoms per unit cell and two different types of atoms. 
For this case n a = 4 and the position vectors of the four basis atoms in the unit cell are 
given by 



di = -(ti + t 2 ) 
d 2 = -(ti - t 2 ) 
d 3 = (ti + t 2 ) 
d 4 = (ti - t 2 ). 



(23) 
(24) 
(25) 
(26) 



The local pseudopotential of the valence electron interacting with a periodic wurtzite 
lattice is given by 



Vl{v) = E 
G 



1 



E^(G)expHG-d j 



exp(iG ■ r) 



(27) 



E — H(G) [exp(-iG • d x ) + expHG • d 4 )] 
G Ua 

+ v 2 (G) [exp(-iG • d 2 ) + exp(-iG • d 3 )]} exp(iG • r), 



(28) 
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since atoms 1 and 4, and 2 and 3 are identical in the wurzite lattice. Rewriting v\ and v 2 in 
terms of Vg and Va from Eq. ( |18|) we obtain 

(29) 



VlW = £ [Vs(G)S s (G) + tV A (G)S A (G)} exp(zG • r) 



G 

where the symmetric (Ss) and anti-symmetric (Sa) structure factors are given by 



S S (G) = —J2 exp(-zG ■ dj), S A (G) = — £ exp(-zG • d, 



n. 



n„ 



(30) 



where Pj = +1 when j = 1,4 and Pj = —1 when j = 2,3. Thus the pseudopotentials of 
wurtzite and zinc-blende crystals differ from each other only in the definition of the structure 
factors. We can carry out the summations in Eq. (|30|) as follows. 

The hexagonal lattice is characterized by three parameters: ao, Co, and uq. Like zinc- 
blende, wurtzite lattice has tetrahedral coordination about each ion, but the orientation 
of the tetrahedron is different from that of zinc-blende. If we assume perfect tetrahedral 
coordination, then 



u = 0.375 



(31) 



£°= /- 

a V 3' 

The direct lattice primitive translation vectors of the simple hexagonal Bravais lattice are 



1 \/3 c 

ai = (l,0,0)a , a 2 = (-, — ,0)a , a 3 = (0, 0, — )a . 

2 2 a 

With this definition, the position vectors of the basis atoms are 



di 
d 2 
d 3 
d 4 



1 r 1 c o 

~1 1 » V 1 + 2 «o)— ao 
4 [ V3 

1 r 1 / Co 

7 1 '^'( 1 + 2m o)— a 

4 [ \/c5 do 

4 [ a/3 a 



(32) 

(33) 
(34) 
(35) 
(36) 



in cartesian coordinate representation. 

In the zinc-blende crystal the nearest neighbor atoms are located half-way along the face 
diagonal of the fee lattice. Consequently, if ao(zb) is the lattice constant of the zinc-blende 
crystal, then the nearest neighbor distance is 



r n (zb) = -^=a (zb). 



(37) 



In the hexagonal crystal, the nearest neighbor atoms are located along the edges of the 
hexagon. Consequently, 



S 



r n (hex) = ao(hex) (38) 
Since r n (zb) = r n (hex), comparing Eqs. fl37D and (^) we obtain the relation 

a (zb) = v^2a (hex). (39) 

and 

V2ir 



G 



ao(hex 



■(G x , G y , G z ). 



(40) 



This definition allows the comparison of G vectors of zinc-blende and wurtzite crystals on 
an equal footing. 

Substituting Eqs. ( j23"l - [26| ) into (|30"D and carrying out some algebraic manipulations we 
obtain 



where 



S S (G) = cos(G ■ ti) cos(G ■ t 2 ), S A (G) = - cos(G ■ ti) sin(G ■ t 2 ) 



tl = ^ (ai + a2) + T = ( I'^I'7! )ao ' 



(41) 



(42) 



t 2 = -u a 3 = (0, 0, v /%/2)a . 



Finally, substituting Eq. (f4Tf ) into ( ¥T|) we obtain 



^^(G) = cos 



y/2 



G nr. G<i, G? 



71 



4 V48 \/6, 



cos(V27rG z y/uo/2) 



(43) 



(44) 



V 4 ^48 \/6, 



sin(v / 27rG, v ^/2). 



(45) 



Sa(G) = — cos 
We can also define G as 

G = (Zbi,mb2,nb 3 ), (46) 
where b 1; b 2 , and b 3 are the reciprocal lattice primitive translational vectors given by 
V2tt 



bi 



a 



:(V2,-^/|,0), 



Substituting Eqs. (pE7|) and fl46|) into (p0|), utilizing the relations 

^ ■ = 27T, b»-a,- = (i ^ j), 
and carrying out some algebraic manipulations, we obtain 



b 2 = — (0,a/-,0), b 3 =— (0,0,^^). (47) 



(48) 
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5*5 (G) = cos 



(I m n 
2tt - + — + - 
\6 6 4 



cos(n7rao) 



S A (G) 



cos 



(I m n N 
2vr - + — + - 
\ 6 6 4, 



sin(n7ruo). 
(49) 



Eqs. Q44]) -(|45[) and (f49|) are equivalent. We find that the above structure factors yield band 
structures in excellent agreement with those of Cohen and Bergstresser for wurtzite crystals 
||39|| . We tested our program with both these structure factors and verified that they give 
identical band structures. 

Previously we used a different coordinate system for the definition of the direct lattice 
vectors of the simple hexagonal Bravais lattice. That definition resulted in structure factors 
different from the ones given above. While these two definitions are equivalent in principle, 
in practice the values of the parameters Vs and Va are intimately linked to the choice of 
the coordinate system and the structure factors determined therefrom. The present choice 
of the coordinate system that gives rise to the structure factors in Eqs. 



D-© and © 

yield correct band structures, in complete agreement with the original band structures of 
Cohen and Bergstresser |39| . 



III. BAND STRUCTURE MODEL 

The virtue of EPM is that it reproduces the bulk band structure to better than 0.1 
eV accuracy [^DJ. Other methods are less accurate, especially for the calculations of the 
band gaps. For example, the calculations based on the density functional theory (DFT) 
typically underestimate the band gaps by about 30-50%. The DFT based methods optimize 
the orbitals of the occupied electronic states only, not those of the unoccupied orbitals. 
This problem remains even when non-local density gradient corrections and other higher- 
order improvements to the DFT methodology are made. Consequently, the band gaps 
determined by DFT are in general grossly in error. Furthermore, the DFT calculations 
are computationally far more expensive compared to EPM. Hence, EPM is suitable for the 
investigation of the electronic structures of semiconductor clusters. 

To apply EPM to the electronic structure calculations of the semiconductor clusters 
we assume that these clusters have the crystal structure of the bulk semiconductor. This 
assumption is justified because we are considering relatively large clusters containing tens to 
hundreds of atoms. Furthermore, the X-ray and TEM experiments on Si, CdS, CdSe, CdTe, 
GaAs, and InP clusters have shown that the bulk lattice structure is preserved even when 
the cluster radii are as small as R = 7 A @,|,|I7]^,[25|] . The reason for the preservation 



of the bulk lattice structure in such small clusters may be due to the presence of ligands 
on the surfaces of these clusters. These ligands are necessary to prevent the clusters from 
coalescing into larger units. These ligands also terminate the dangling bonds on the surfaces 
of these clusters, and thus inhibit structural reconstruction. For this reason, even small 
clusters seem to possess bulk lattice structure. The major effect of size on cluster structures 
seems to be small contractions of the bonds relative to their bulk bond lengths. 

The band structure calculations are carried out for these clusters in almost the same 
way as we had done previously for CdS [30|. In bulk semiconductors the allowed wave 



vectors k of the electronic states are continuous. On the other hand, only discrete k- 
states are allowed in clusters. If we model the cluster as a rectangular box with dimensions 
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L x ,L y , and L z , then a reasonable approximation is to use the bulk pseudopotential V p (r) 
inside the cluster and terminate this potential at the surfaces of the cluster by an infinite 
potential. The wave vectors of the lowest allowed states are then given by the quantization 
condition sin(k x L x ) sm(k y L y ) sin(k z L z ) = 0, whose solution is 

k = 7rfe^V (50) 

where n x , n y , and n z are the quantum numbers of a particle in a box with infinite potentials 
at the boundaries. The surface ligands act a potential barriers to the valence electrons. 
Consequently, for low energy excitations under consideration, the assumption of infinite 
potentials at the boundaries is a good approximation. The energy levels at these allowed 
k-states constitute the band structure of a rectangular box. 

Similarly, if we model the cluster as a spherical object of radius R, the energy levels of 
the valence electrons will be quantized because of the spherical boundary. The wave vectors 
of the lowest allowed states are given by j (k n R) = 0, whose solution is k n = tvk/R p0| . 
Since k n is along the radial direction, we project it onto each of the cartesian axes with 
equal magnitude to obtain cartesian components of k. This procedure yields 

7T 

(n x ,n y ,n z ). (51) 



V3R 



for spherical clusters. Sometimes we also model spherical clusters as cubic boxes with L = 2R 
as the sidelength of the cube. The energy levels at these allowed k-states constitute the band 
structure of a spherical cluster. 

The exciton energy of a cluster of radius R is given by [[30 



1 786 

E x = E g - — - 0.248£^, (52) 

where E g is the band gap, e is the dielectric constant, E Ry = fxe 4 /2e 2 h 2 is the effective 
Rydberg energy of the exciton, and fi is the reduced mass of the electron-hole pair. 

EPM has been shown to reproduce the band gaps to within 0.1 eV for both bulk materials 
and clusters J30J. Specifically, our predicted indirect band gap of 0.43 eV for a 10 A radius 



silicon cluster was found to be in good agreement with recent experimental value of 0.5 eV 
obtained by Louis Brus and co-workers ||41||. Likewise, EPM yielded excellent results for CdS 



clusters in comparison with experiment |HJ . Full details on the computational methodology 



are given in Refs. . The utility of the proposed method for the calculations of the spectral 



shifts is also documented in the previous publications pi 



IV. COMPUTATIONAL DETAILS 

The binary semiconductors being considered in this study all consist of either zinc- 
blende or wurtzite lattices. Consequently, the Hamiltonian matrices are of complex valued 
Hermitian matrices. We diagonalized these matrices using the EISPACK routines CH, 
HTRIDI, and TQLRAT. Typically, we use 283 plane waves to converge the energies to 
better than 0.01 eV accuracy. 
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Table I lists the parameters of the local pseudopotentials of zinc-blende and wurtzite CdS 
and CdSe clusters, while Tables II and III list the corresponding parameters for the non- 
local pseudopotentials of zinc-blende CdTe, A1P, GaP, GaAs, and InP clusters. Spin-orbit 
interaction was included in the calculations of wurtzite CdS, wurtzite CdSe, zinc-blende 
CdTe, and InP clusters. We have verified that the spin-orbit effects are small in A1P, GaP, 
and GaAs crystals. Consequently, our calculations on these clusters omitted the spin-orbit 
interaction in the Hamiltonian. 

The s orbitals don't exhibit spin-orbit interaction since they have I = orbital angular 
momentum. The innercore and d electron states are eliminated in EPM through the use of 
the pseudopotentials. Consequently, only the outermost p orbitals are affected by the spin- 
orbit interaction. Hence, we used the radial Hartree-Fock-Slater orbitals 4p for Cadmium, 
5p for Indium, 3p for Sulphur, 4p for Selenium, 5p for Tellurium, and 3p for Phosphorus 
in our calculations employing spin-orbit interaction in the Hamiltonian. Other calculations 
reported in the literature seem to include innercore p orbital wave functions in the Hamilto- 
nian, instead of the outermost p orbitals we employ. However, we have verified that in each 
case we obtain better band structures with the outermost p orbital wave functions than with 
the innercore p orbitals. 

We use the standard notation to represent the high symmetry points in the Brillouin 
zone P0| . For zinc-blende crystals 

X = (1, 0, 0), W = (1, 0.5, 0), L = (0.5, 0.5, 0.5), F = (0, 0, 0), U = K = (0, 0.75, 0.75) (53) 

in units of l-nja^. For wurtzite crystals, the high symmetry points are defined as 

r = (0, 0, 0), K = (V2/3, V2/3, 0), M = (V2/3, 0, 0), 

A — (0, 0, a )/\/2co) = (0,0,^/4), 

H = \/2A Vv^co) = (^2/3, y/2/3, v/3/4), 

L = {y/2/3, 0, a /V2co) = (y/2/3, 0, v/3/4) (54) 

in units of ^/2tc/clo- We map bands along these symmetry points to obtain the complete 
band structure. 

V. RESULTS AND DISCUSSION 

A. CdS and CdSe Clusters 

The CdS and CdSe crystals exist in both zinc-blende and wurtzite form, with the wurtzite 
form being the ground state structure. In principle it is possible to synthesize clusters also in 
both these lattice forms. However, the CdS clusters seem to prefer zinc-blende over wurtzite, 
while CdSe clusters seem to prefer wurtzite over zinc-blende structure. The possibility of 
being able to synthesize both these isomeric forms of CdS and CdSe clusters, provides an 
opportunity to investigate the effect of lattice structure on the exciton energies as a function 
of cluster size. 

We use the lattice constant a = 5.82 A for zinc-blende and = 4.14 A for wurzite 
CdS crystals. Our calculated band gaps are 2.44 eV and 2.52 eV for bulk zinc-blende and 
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wurtzite CdS, respectively. The experimental result is 2.50 eV [0 for both these structures. 
We use the lattice constant a = 6.08 A for zinc-blende and a = 4.30 A for wurzite CdSe 
crystals. The calculated band gaps for CdSe are 1.92 eV and 1.79 for bulk zinc-blende and 
wurtzite CdSe, respectively. The experimental values for bulk wurtzite CdSe are 1.75 and 
1.83 eV at 80 and 293 K, respectively. We take the average of these two values, 1.79 eV, 
as the experimental value for both zinc-blende and wurtzite CdSe crystals. According to 
Cohen and Bergstresser the splitting in energy of the fundamental gap without spin- 
orbit coupling (Tg — > T^) is E a + |A, where E a is the value of the band gap with spin-orbit 
coupling (Tg — > T^) and A is the value of the spin-orbit splitting. The spin-orbit coupling 
parameter [i was chosen in order to satisfy this condition to better than 0.001 eV accuracy. 
For CdS crystal, A = 0.062 eV at 77K |43| and we determined that /i = 0.00008 is optimal. 
For CdSe crystal, A = 0.42 eV at 77 K ||3j and we determined that \x = 0.00041 is optimal. 
Figure 1(a) show the band structure of a wurtzite CdS crystal, including the spin-orbit 
coupling. The discrete energy levels of a R = 15.0 A zinc-blende CdS cluster, modelled as 
a sphere of radius R, are given in Fig. 1(b). Figure 1(c) show the band structure of a zinc- 
blende CdS zinc-blende crystal, compared with experimental data |[44| . Figure 2(a) show 
the band structure of a wurtzite CdSe crystal, including the spin-orbit coupling, and 2(b) is 
the band structure of a zinc-blende CdSe crystal, compared with the available experimental 
H|. The study of the dependence of the excitonic energies on the cluster size was 

5.5, m e = 0.19, m h = 0.80 @ for both 
0.13, m h = 0.45 OT for both 



data 



carried out using the bulk crystal parameters e 
wurtzite and zinc-blende CdS system, and e = 



30 



10.0, m e = 0.13, m h = 0.45 
wurtzite and zinc-blende CdSe system. The exciton energies are reported in Tables IV- VII. 
Since for a cluster of radius R = 100 A the fundamental gap is 2.46 eV for zinc-blende, 
and 2.55 eV for wurtzite, we shifted our conduction energy levels by +0.04 eV and -0.05 
eV, respectively, in order to obtain the experimental bulk band gap of 2.50 eV. We carried 
out a similar correction of for CdSe clusters, shifting the conduction bands by -0.15 eV for 
zinc-blende, and -0.025 eV for wurtzite, so that we obtain the bulk band gap of 1.79 eV 
for a cluster of R = 100 A. In Fig. 3 we compare the calculated exciton energies of both 
wurtzite and zinc-blende CdS clusters with the experimental data [p§l , f45| -48l. This figure 



clearly shows that the zinc-blende and wurtzite clusters exhibit different spectral shifts and 
the available experimental data follow one or the other of these trends. Based on these 
calculations, we can assign the zinc-blende structure to the clusters synthesized by Wang 
and Herron [47[] and wurtzite structure to the smaller clusters synthesized by Weller and 
coworkers IktHI. 



1. The Effect of Dielectric Constant and Spin- Orbit Coupling 

We have investigated the dependence of the exciton energies on dielectric constant (e) 
by carrying out calculations on CdSe clusters using two different values of e: e(0) and e(oo). 
Figure 4(a) compares the calculated exciton energies with the available experimental data 
P9H54]. The dashed line represents the calculations on wurtzite CdSe clusters with e(0) = 
10.0, while the solid line represents the corresponding calculation with e(oo) = 6.25. Accord- 
ing to these results the exciton energies calculated with e(oo) are closer to the experimental 
data. Figure 4(b) represents identical calculations on zinc-blende structure, with m e = 0.11 
and rrih = 0.44 |41|, compared with the experiments pl|-|5^||. In this case the exciton energies 
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seem to be insensitive to a reasonable choice of the dielectric constant. 

The bulk wurtzite CclSe semiconductor exhibits three distinct transitions, labelled A, 
B, and C, arising from the splitting of the valence band due to spin-orbit interaction. If 
we assume that similar transitions will exist in clusters too, then we obtain the transitions 
shown in Fig. 4(c). Since the bulk transition A is at 1.79 eV |}4"3"|| , and that for a cluster of R 
= 100 A is at 1.67 eV (using 233 plane waves), a correction of +0.12 eV was applied. The 
calculated transitions B and C were also similarly corrected: B bulk transition = 1.81 eV 
cluster (R = 100 A) = 1.74 eV, so the correction = 0.07 eV; C bulk transition = 2.22 eV 



4*51 , cluster (R = 100 A) = 2.14 eV, so the correction = 0.08 eV. From Fig. 4(c) we observe 
that if only the lowest exciton energies are considered, then the experimental data show 
significant scatter compared to the calculated results. However, when we consider all three 
transitions (A, B, and C), then all the experimental data can be neatly accounted for as 
belonging to one of these three transitions. This finding indicates the possibility that some 
of the experimental data reported in the literature correspond to higher-energy transitions 
(B and C) rather than to the lowest energy transition A. 

The spectroscopy of nanoscale CdSe clusters have attracted considerable attention from 
several research groups |49|-">4[] . Some of these groups have succeeded in synthesizing CdSe 
cluster samples that have very narrow size distributions |5l|JiS3|j54]] . Furthermore, they have 
made very careful measurements of the exciton energies. For these reasons, direct comparison 
between our calculations and the data from these groups provides a benchmark test of 
the reliability of our computational method in yielding accurate exciton energies. From 
the results presented above it is clear that the agreement between our calculations and 
experiment is excellent. The experimental data lie much closer to our calculated values than 
to those of the EMM, without any exception. The experimental measurements on cluster 
samples with very narrow size distribution are in better agreement with our calculations 
than the corresponding data on samples with broad size distribution. Finally, we are able to 
identify that some measured exciton energies probably correspond to higher energy B and 
C transitions than to the lowest energy A transition. 



B. CdTe Clusters 

Since both Cadmium and Tellerium are large atoms, the CdTe clusters and crystals exist 
only in zinc-blende form. Since both the cation and anion are heavy atoms, we expect 
significant spin-orbit coupling in CdTe. Consequently, we have carried out band structure 
calculations on zinc-blende CdTe clusters using non-local EPM with the effects of spin-orbit 
coupling included in the Hamiltonian J5f*J]. The parameters used in these calculations are 
given in Tables II and III. The spin-orbit coupling parameter /i was fit to the experimental 
band gap of 1.56 eV at 300 K g|. 



Figure 5(a) presents the band structure for the zinc-blende CdTe crystal, while Fig. 5(b) 
compares the calculated direct exciton energies with the available experimental data |6li|62| . 
In these calculations we used e = 10.2, m e = 0.09, and m h = 0.72 JJ3[. The exciton energies of 



spherical clusters of different radii are reported in Table VIII. The computational procedure 
employed is identical to that described above. At present, experimental data on CdTe 
clusters are few. However, we hope that our calculations will stimulate further experimental 
efforts on these clusters. 
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C. A1P Clusters 



The binary semiconductor A1P is isoelectronic to Si2 and they both crystallize in an 
identical lattice structure. Also, since Al and P are the nearest neighbors of Si in the 
periodic table, it is reasonable to expect that crystalline A1P will have electronic properties 
similar to that of bulk Si. Indeed that turns out to be true: both Si and A1P are indirect 
gap semiconductors with similar band structures. Likewise, since both these semiconductors 
are made of light elements, we expect spin-orbit coupling to be negligible in these crystals. 
We have reported on our investigation of the spectral shifts of silicon clusters before. Now 
we present corresponding calculations on A1P clusters. 

The previous calculated band structure of A1P was obtained using the LCAO-MO method 
and the corresponding band structure is not available from EPM calculations. Consequently, 
we have fitted the non-local parameters of A1P to reproduce the band structure calculated 
using the LCAO-MO method. Our initial guess for these parameters utilized the parameters 
of bulk Si fSOfl . In the end the optimal parameters gave 3.63 eV for the direct and 2.41 eV for 
the indirect gaps. These calculated values are in complete agreement with the corresponding 
experimental band gaps at 300 K [p^j . Tables II and III give the parameters, while Fig. 6 
gives the bulk band structure and the calculated exciton energies of A1P clusters. These 
calculations utilized e = 9.8, m e = 0.21, and m h = 0.94 The exciton energies are 

reported in Table IX. Nanoscale A1P clusters have not been synthesized in the laboratory 
so far. However, recent experimental efforts on Si clusters point to the possibility of similar 
interest on A1P clusters. 



D. GaP and GaAs clusters 

Figures 7 and 8 give the band structures and the exciton energies of GaP and GaAs 
spherical clusters as a function of their radius. Tables II and III show the parameters that 
were used in these calculations |5"5| . Figures 7 (a) and 8 (a) show the band structures of 
zinc-blende GaP and GaAs crystals obtained using the non-local empirical pseudopotential 
(solid line) and the local empirical pseudopotential (dashed line). Comparison with the 
experimental data of GaP [^,0 and GaAs |]66|-|69|| crystals clearly shows that the non-local 
approximation yields to band structures that are better in agreement with the experimental 
data than those one obtained with the local approximation. Since GaP is an indirect gap 
semiconductor, we present spectral shifts of both the direct and indirect transitions. The 
results obtained from the use of both local and non-local EPM are shown in these figures. 
The calculated band gaps of GaP, obtained using 283 plane waves and the local EPM 
method, are 2.79 eV for the direct transition and 2.15 eV for the indirect transition. Since 
the corresponding experimental band gaps are 2.78 eV for the direct transition and 2.27 eV 
for the indirect transition [|55| , we shifted the calculated band gaps by -0.01 eV, and +0.12, 
respectively. The non-local EPM gives 2.88 eV and 2.17 eV, respectively, for the direct 
and the indirect band gaps of GaP. Consequently, the corresponding shifts are -0.10 eV and 
+0.10 eV. 

For GaAs, the experimental direct gap is 1.47 eV, which is an average of the experimental 
band gaps at and 300 K |^3| . The calculated band gaps, obtained using 283 plane waves, are 



1.50 eV with the local calculation and 1.52 eV with the non-local calculation. Consequently, 



15 



the corrections are -0.03 eV and -0.05 eV, respectively. We used the bulk crystal parameters 



e = 9.1, m e = 0.10, m h = 0.86 || for GaP clusters, and e = 9.1, m e = 0.07, m h = 0.86 gO 
for GaAs clusters. 

The exciton energies thus calculated are presented in Tables X-XII and Figs. 7 and 
8. From these figures it is clear that over a large range of cluster sizes the local EPM is 
able to reproduce spectral shifts as accurately as non-local EPM for both GaP and GaAs. 
However, at small cluster sizes the non-local correction on the spectral shifts is significant. At 
large cluster sizes the absorption spectrum shifts to higher energies with decreasing cluster 
size. This blueshift is expected due to confinement of the electron-hole pair in the cluster. 
However, at small cluster sizes the absorption spectrum of GaAs clusters shifts to lower 
energies with decreasing cluster size; a trend opposite to that observed for large clusters. 
In the case of GaP, which is an indirect gap semiconductor, the lowest energy transition 
exhibits the expected blueshift at both large and small cluster sizes. But this transition is 
not observable because it is forbidden. The origin of the absorption spectrum, corresponding 
to the observable direct transition, shifts to lower energies with decreasing cluster size at 
small cluster sizes. Both local and non-local pseudopotentials exhibit the same qualitative 
behavior. The main difference is that the non-local EPM predicts less redshift in small 
clusters than the local EPM. 

We can explain the calculated trends in the following way. At large cluster sizes the 
electron and hole are both confined in a spherical well. This quantum confinement increases 
the band gap with decreasing cluster size and it is the dominant effect in this size regime. 
In these large clusters, the negatively charged electron and the positively charged hole are 
spatially separated and hence the coulomb attraction between them is negligible. However, 
in small clusters the coulomb attraction energy between the electron-hole pair cannot be 
neglected. While the band gap still increases with decreasing cluster size, in small clusters 
this increase is sufficiently overcome by the coulomb energy that the spectra shift to lower 
energies. Consequently, in this small cluster size regime the absorption spectra of clusters 
may exhibit redshift instead of the blueshift. 

At present reliable experimental data are not present for the spectral shifts of these im- 
portant III-V semiconductor clusters, partly because of considerable experimental difficulties 
that arise in trying to synthesize III-V semiconductor clusters with narrow size distribution. 
However, based on our experience with CdS and CdSe clusters, and InP clusters (see below), 
we expect our calculated spectral shifts of A1P, GaP, and GaAs clusters to be good estimates 
of the expected spectral shifts. 



E. InP clusters 



Recently, Nozik and co-workers have succeeded in synthesizing InP clusters and for the 
first time showed an exciton transition in the spectrum of a III-V quantum dot |7D |. The 



present work compares their experimental results with our theoretical predictions. We used 
the same procedure as before: First, we found the parameters of the EPM Hamiltonian and 
then calculated the excitonic energies of InP clusters. The band structure calculations on 
zinc-blende InP crystal were carried out using the non-local EPM that includes the effects 
of spin-orbit coupling. The parameters we employ are slightly different from those used by 
Chelikowsky, but otherwise the two calculations give nearly identical band structures. In 
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particular, the coupling parameter // = 0.00023 was chosen in order to reproduce the fun- 
damental splitting = 1.50 eV, as obtained by Chelikowsky ||60|| . Complete set of parameters 
used in our calculations are tabulated in Tables II and III. 

Figure 9(a) displays the band structure of bulk InP and Fig. 9(b) presents the exciton 
energies of the clusters as a function of their radius. These calculations have employed the 
parameters e = 9.61, m e = 0.077, and nih = 0.58 |63[| . We shifted all the conduction bands 
by -0.207 eV so that the band gap of an R = 100 A cluster is equal to the experimental 
bulk value of 1.35 eV at 300 K. The exciton energies are reported in Table XIII. 

The type (a) sample in the experiments of Nozik and co-workers have a radius of about 
20 A and exhibits an exciton transition at about 1.7 eV ||70| . This paper quotes 0.35 
eV blueshift for the onset of absorption for the colloidal sample (a) whose mean diamter 
is 26.1 A. However, 0.35 eV is not the blueshift for a particle with a diameter of 26.1 
A because colloid (a) has a broad size distribution with a standard deviation of 7.5 A (see 
Fig. 3 of Ref. J7DJ ). Thus the onset of absorption really corresponds to a size of about 
40 A (private communication). We use this value for the particle size while comparing 
experimental blueshift with our calculations. The difference between our calculated exciton 
energy and the experimental value is approximately 0.17 eV. The type (c) sample has a 
radius of 13 A and exciton transition at 2.25 eV, in perfect agreement with our calculation. 
From this comparison it is clear that the overall trend of the experimental spectral shifts is 
in close agreement with our calculations. Additional reliable data on spectral shifts over a 
range of cluster sizes will certainly be useful in establishing the accuracy of our calculations. 

Our band structure calculations on all IIIA-VB semiconductor clusters have shown that 
local empirical pseudopotentials are reasonably accurate compared to the non-local empirical 
pseudopotentials for the calculations of the spectral shifts of these clusters. The non-local 
corrections on the spectral shifts are most important in small cluster sizes. In addition, 
our calculations have shown that, while quantum confinement energy is the dominant factor 
affecting spectral shifts in large clusters, the Coulomb interaction between the electron and 
hole has significant effect in small clusters. The attractive Coulomb interaction is sufficiently 
strong in small clusters that it overcomes the confinement energy of the electron-hole pair 
and gives rise to redshift, instead of the blueshift, of the electronic absorption spectrum. 

VI. SUMMARY 

In summary, we calculated the spectral shifts of CdS, CdSe, CdTe, A1P, GaP, GaAs, and 
InP semiconductor clusters using the most accurate available empirical pseudopotentials. 
These semiconductors cover a wide range of bond polarities and band structures. These 
binary clusters also represent a series in which one ion is held constant while the other ion 
is varied along a column of the periodic table. Furthermore, the pseudopotentials employed 
in these calculations incorporate the effects of non-locality and spin-orbit coupling when- 
ever they are important. For many of these semiconductors, we had to first determine the 
pseudopotential parameters based on the latest experimental and theoretical data. Conse- 
quently, we report complete set of pseudopotential parameters employed in our calculations. 
Furthermore, we also give full details of our band structure calculations employing pseudopo- 
tentials that incorporate both non-local and spin-orbit effects. At present these calculations 
represent the most sophisticated calculations on the spectral shifts of semiconductor clusters. 
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Previously we had shown that a simple local pseudopotential yields spectral shifts of 
zinc-blende CdS clusters in excellent agreement with experiment over a range of cluster 
sizes. Now we extend these calculations to include the effects of spin-orbit coupling in 
the pseudopotential. As before, we find that for zinc-blende CdS, wurtzite CdS, wurtzite 
CdSe, and zinc-blende InP clusters the spectral shifts calculated using our band structure 
model are in excellent agreement with experiment. The shapes and crystal structures of 
the unit cell have significant effect on the exciton energies. The small clusters in particular 
are sensitive to whether their crystal structure is zinc-blende or wurtzite. In the case of 
small CdS clusters where the experimental data on lattice structure are either ambiguous or 
unavailable, our calculations are able to assign the structure unambiguously. While reliable 
experimental data are not yet available at present for zinc-blende CdSe clusters, we predict 
that the spectral shifts in these clusters will be nearly the same as those of wurtzite clusters. 

In the absence of experimental data on CdTe, A1P, GaP, and GaAs clusters, our calcula- 
tions provide reasonable estimates of the expected spectral shifts and trends in these clusters 
as a function of cluster size. The very little difference between the local and the non-local 
calculations confirms the validity of the local pseudopotential method for the calculations 
of the exciton energies. The main effect of the non-local pseudopotentials seems to be to 
reduce the magnitude of the redshift seen in the spectral shifts of some small III-V semi- 
conductor clusters. Similarly, the spin-orbit interaction does not change the lowest exciton 
energies. The main effect of the spin-orbit interaction is to split the valence band into sub- 
bands, thus giving rise to new transitions in the spectra. Some experimental data seem to 
be in better agreement with these higher energy transitions, thus suggesting the possibility 
that transitions originating from spin-orbit split valence bands are being observed even in 
these nanoscale clusters. This also conceivably explains the reasons for considerable scatter 
in the experimentally determined exciton energies for a given cluster size. This study, to- 
gether with our previous investigation, provides evidence that the empirical pseudopotential 
method yields unique insights into the quantum confinement effects and is a powerful tool 
for calculating the spectral shifts of semiconductor clusters. 
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TABLES 



TABLE I. Reciprocal lattice vectors and form factors (in a.u.) for CdS and CdSe crystals. 
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TABLE II. The Pseudopotential parameters of the semiconductors. The parameters Vs, Va, 
ao, and A2 are in a.u., ao is in A, and 0q and [i are dimensionless. 



Compound 


V S (V3) 


V S (V8) 


v s (Vn) 


v A (Vs) 


Va(V4) 


v A (Vn) 


a (A) 


Local parameters 
















CdTe 


-0.11000 


0.00000 


0.03100 


0.03000 


0.02500 


0.01250 


6.48 


A 1 T ~i 

A1P 


-0.12250 


r\ r\ -\ r\c\r\ 

0.01000 


0.02000 


0.05000 


f\ A'T AAA 

0.03000 


l~i AA AA r 

0.00005 


5.46 


GaP 


-0.11500 


0.01000 


0.02900 


0.05000 


0.03500 


0.01300 


5.45 


GaAs 


-0.10700 


0.00700 


0.03900 


0.02800 


0.01900 


0.00100 


5.65 


InP 


-0.11750 


0.00000 


0.02650 


0.04000 


0.03000 


0.01500 


5.86 


Non-local parameters 




















Cation 






Anion 




Spin Orbit 




OlQ 


00 


A 2 




00 

/ U 


A 2 




CdS 














0.00008 


CdSe 














0.00041 


CdTe 


0.000 


0.400 


0.000 


0.000 


0.400 


1.000 


0.00061 


A1P 


0.190 


0.300 


0.350 


0.160 


0.030 


0.180 




GaP 


0.000 


0.300 


0.200 


0.160 


0.050 


0.225 




GaAs 


0.000 


0.000 


0.063 


0.000 


0.000 


0.313 




InP 


0.000 


0.250 


0.275 


0.150 


0.050 


0.175 


0.00023 



TABLE III. The Pseudopotential parameters Rq and R2 in A. 



Compound 



Cation 



Anion 



CdTe 
A1P 
GaP 

GaAs 
InP 



Ro 
1.37 
1.06 
1.27 
1.27 
1.27 



R2 
1.41 
1.19 
1.19 
1.44 
1.29 



Ro 
1.06 
1.06 
1.06 
1.06 
1.06 



R2 
1.41 
1.19 
1.19 
1.44 
1.29 
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TABLE IV. Band gaps and exciton energies of zinc-blende CdS spherical clusters obtained 
using 283 plane waves. R is the radius of the cluster, k is the angular momentum wave vector 
of the lowest exciton state, E g is the band gap, V c is the Coulomb energy, and E x is the exciton 
energy. 



R{A) fc(27r/a ) = ap/VUR E g (eV) V c (eV) E x (eV) 

5.0 0.3359 4.69 -0.93 3.77 

10.0 0.1680 3.69 -0.46 3.24 

15.0 0.1120 3.14 -0.31 2.85 

30.0 0.0560 2.65 -0.15 2.51 

50.0 0.0336 2.52 -0.09 2.44 

99.0 0.0170 2.46 -0.04 2.43 



TABLE V. Band gaps and exciton energies of wurtzite CdS spherical clusters obtained using 
233 plane waves and neglecting the spin-orbit coupling. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E g is the band gap, V c is the Coulomb 
energy, and E x is the exciton energy. 



R{A) k(V2ir/a ) = oq/VGR E g (eV) V c (eV) E x (eV) 

5.0 0.337700 5.66 -0.93 4.66 

10.0 0.1689 3.92 -0.46 3.39 

15.0 0.1126 3.27 -0.31 2.89 

30.0 0.0563 2.74 -0.15 2.52 

50.0 0.0338 2.60 -0.09 2.45 

99.0 0.0171 2.54 -0.04 2.43 



TABLE VI. Band gaps and exciton energies of zinc-blende CdSe spherical clusters obtained 
using 283 plane waves. R is the radius of the cluster, k is the angular momentum wave vector 
of the lowest exciton state, E g is the band gap, V c is the Coulomb energy, and E x is the exciton 
energy. 



R(A) fc(27r/a ) = a /Vl2R E g (eV) V c (eV) E x (eV) 

5.0 0.351000 4.18 -0.51 3.51 

10.0 0.1755 3.18 -0.25 2.77 

15.0 0.1170 2.64 -0.17 2.31 

30.0 0.0585 2.14 -0.08 1.90 

50.0 0.0351 2.00 -0.05 1.79 

99.0 0.0177 1.94 -0.02 1.76 
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TABLE VII. Band gaps and exciton energies of wurtzite CdSe spherical clusters obtained using 
233 plane waves and neglecting the spin-orbit coupling. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E g is the band gap, V c is the Coulomb 
energy, and E x is the exciton energy. 



R{A) k(V2ir/a ) = a /V6R E g (eV) V c (eV) E x (eV) 

5.0 0.3510 5.00 -0.51 4.45 

10.0 0.1755 3.28 -0.25 3.00 

15.0 0.1170 2.60 -0.17 2.40 

30.0 0.0585 2.04 -0.08 1.92 

50.0 0.0351 1.88 -0.05 1.80 

99.0 0.0177 1.81 -0.02 1.76 



TABLE VIII. Band gaps and exciton energies of zinc-blende CdTe spherical clusters obtained 
using 137 plane waves and including the spin-orbit coupling. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E g is the band gap, V c is the Coulomb 
energy, and E x is the exciton energy. 



R{A) k(27r/a ) = ap/VUR E g (eV) V c (eV) E x (eV) 

5.0 0.3741 3.40 -0.50 2.89 

10.0 0.1871 3.00 -0.25 2.75 

15.0 0.1247 2.53 -0.16 2.36 

30.0 0.0624 1.94 -0.08 1.86 

50.0 0.0374 1.75 -0.05 1.69 

99.0 0.0189 1.65 -0.02 1.62 



TABLE IX. Band gaps and exciton energies of zinc-blende A1P spherical clusters obtained 
using 283 plane waves. R is the radius of the cluster, k is the angular momentum wave vector of 
the lowest exciton state, E x & is the direct exciton energy, and E X i is the indirect exciton energy. 



R{A) k(27r/a ) = a /Vl2R E xd (eV) E xi (eV) 

5.0 " 0.3154 3.91 " 3.83 

10.0 0.1577 4.02 2.93 

15.0 0.1051 3.97 2.68 

30.0 0.0526 3.73 2.46 

50.0 0.0315 3.63 2.41 

99.0 0.0159 3.58 2.39 
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TABLE X. Band gaps and exciton energies of zinc-blende GaP spherical clusters obtained 
using 283 plane waves and the local pseudopotential. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E xd is the exciton energy for the direct 
transitions, and E x i is the exciton energy for the indirect transitions. 



R{A) k(27r/a ) = ap/VUR E xd (eV) E xi (eV) 

6.0 0.2624 3.09 3.15 

12.0 0.1312 3.43 2.58 

18.0 0.0875 3.25 2.43 

30.0 0.0525 2.99 2.31 

51.0 0.0309 2.84 2.26 

99.0 0.0159 2.79 2.25 



TABLE XI. Band gaps and exciton energies of zinc-blende GaP spherical clusters obtained 
using 283 plane waves and the non-local pseudopotential. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E xd is the exciton energy for the direct 
transitions, and E x i is the exciton energy for the indirect transitions. 



R{A) fc(2vr/a ) = a /VT2R E xd (eV) E xl (eV) 

6.0 0.2624 3.50 3.36 

12.0 0.1312 3.57 2.67 

18.0 0.0875 3.33 2.49 

30.0 0.0525 3.05 2.36 

52.0 0.0303 2.89 2.31 

98.0 0.0152 2.84 2.30 



TABLE XII. Band gaps and exciton energies of zinc-blende GaAs spherical clusters obtained 
using 283 plane waves. R is the radius of the cluster, k is the angular momentum wave vector of 
the lowest exciton state, E x \ is the exciton energy for the direct transitions in a local calculation, 
and E xn i is the direct exciton energy in a non-local calculation. 



R(A) fc(27r/a ) = a /Vl2R E xl (eV) E xnl (eV) 

6.0 0.2718 2.29 2.42 

12.0 0.1359 2.45 2.40 

18.0 0.0906 2.14 2.11 

30.0 0.0544 1.76 1.75 

51.0 0.0320 1.54 1.54 

99.0 0.0165 1.44 1.44 
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TABLE XIII. Band gaps and exciton energies of zinc-blende InP spherical clusters obtained 
using 137 plane waves and including the spin-orbit coupling. R is the radius of the cluster, k is the 
angular momentum wave vector of the lowest exciton state, E g is the band gap, V c is the Coulomb 
energy, and E x is the exciton energy. 



R{A) fc(27r/a ) = ap/VUR E g (eV) V c (eV) E x (eV) 

5.0 0.3384 3.12 -0.53 2.37 

10.0 0.1692 2.89 -0.26 2.42 

15.0 0.1128 2.50 -0.17 2.12 

30.0 0.0564 1.90 -0.08 1.60 

50.0 0.0338 1.67 -0.05 1.41 

99.0 0.0171 1.55 -0.02 1.32 
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FIGURES 



FIG. 1. a) Band structure of a wurtzite CdS crystal implementing the spin-orbit interactions 
and 135 plane waves, b) Allowed electronic levels of wurtzite CdS spherical clusters ( radius = 15 
A) obtained using 233 plane waves, c) Band structure of a zinc-blende CdS crystal obtained using 
283 plane waves. Experimental data are superimposed on the band structure for comparison [44]. 

FIG. 2. a) Band structure of a wurtzite CdSe crystal implementing the spin-orbit interactions 
and 135 plane waves, b) Band structure of a zinc-blende CdSe crystal obtained using 283 plane 
waves. Experimental data are superimposed on the band structure for comparison [44]. 

FIG. 3. Direct exciton energy of wurtzite (upper line), and zinc-blende (lower line) CdS 
spherical clusters compared with the experiments. Filled diamonds are used for the zinc-blende 
data points while filled circles represent the wurtzite clusters. 233 plane waves were used for the 
wurtzite and 283 plane waves for the zinc-blende structures. 

FIG. 4. a) Direct exciton energies of wurtzite CdSe clusters obtained using two different dielec- 
tric constants: e(0) (dashed line), and e(oo) (solid line). Comparison is made with the available 
experimental data: stars [49]; plus [50]; squares [51]; up triangles (light-hole data from [52]); down 
triangles (split-off- hole data from [52]); diamonds [53]; circles (SAXS data from [54]); right tri- 
angles (TEM data from [54]). b) Direct exciton energies of zinc-blende CdSe clusters obtained 
using two different dielectric constant: e(0) (dashed line), and e(oo) (solid line). The comparison 
is made with the available experimental data: up triangles (SAXS data from [54]); down triangles 
(TEM data from [54]); circles [55]; squares [56]; star [57]; plus [58]; diamonds [59]. c) Calculated 
A (lower line), B (intermediate line), C (upper line) transition energies, as discussed in the text. 
The experimental data are represented as in (a). 

FIG. 5. a) Band structure of zinc-blende CdTe crystals obtained using 137 plane waves and 
the spin-orbit interactions, b) Direct exciton energies of zinc-blende CdTe clusters compared with 
the experiments: circles [61]; diamonds [62]. 

FIG. 6. a) Band structure of a zinc-blende A1P crystal obtained using 283 plane waves, b) 
Direct exciton energies of zinc-blende A1P clusters obtained using 283 plane waves, c) Indirect 
exciton energies of zinc-blende A1P clusters obtained using 283 plane waves. 
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FIG. 7. a) Band structure of a zinc-blende GaP crystal obtained using 283 plane waves.. It is 
shown a comparison between the non-local (solid line) and the local (dashed line) empirical approx- 
imation. Experimental data are superimposed on the band structure: diamonds [66]; plus symbols 
[67]. The non-local pseudopotential calculations are seen to agree better with the experimental 
data than the local pseudopotential calculations, b) Direct exciton energies of zinc-blende GaP 
clusters. The open circles are obtained using local pseudopotentials and filled circles are obtained 
using non-local pseudopotentials. c) Indirect exciton energies of zinc-blende GaP clusters. The 
open circles are obtained using local pseudopotentials and filled circles are obtained using non-local 
pseudopotentials . 

FIG. 8. a) Band Structure of a zinc-blende GaAs crystal. It is made a comparison between the 
non-local (solid line), and the local (dashed line) empirical approximation. Experimental data are 
superimposed on the band structure: diamonds [66]; plus symbols [67]; x symbols [68]; squares [69]. 
The non-local pseudopotential calculations agree better with the experimental data than the local 
calculations, b) Direct exciton energies of GaAs spherical clusters. The open circles are obtained 
using local pseudopotentials and filled circles are obtained using non-local pseudopotentials. The 
number of plane waves used is 283. 

FIG. 9. a) Band structure of a zinc-blende InP crystal obtained using 137 plane waves and 
spin-orbit interactions, b) Direct exciton energies of zinc-blende InP clusters compared with ex- 
periments: circles [70]. 
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